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■ 

f^ . We interpret measurements of the Reynolds number dependence of the 

Cn , torque in Taylor-Couette flow by Lewis and Swinney [Phys. Rev. E 59, 

C^l 5457 (1999)] and of the pressure drop in pipe flow by Smits and Zagarola, 

r^ ' [Phys. Fluids 10, 1045 (1998)] within the scaling theory of Grossmann and 

^^ ■ Lohse [J. Fluid Mech. 407, 27 (2000)], developed in the context of thermal 

CN . convection. The main idea is to split the energy dissipation into contributions 

^ from a boundary layer and the turbulent bulk. This ansatz can account for 

Q I the observed scaling in both cases if it is assumed that the internal wind 

Q^ ' velocity U^ introduced through the rotational or pressure forcing is related 

PJ ■ to the the external (imposed) velocity U, by Uw/U oc Re^ with ^ = —0.051 

i-H . and ^ = —0.041 for the Taylor-Couette {U inner cylinder velocity) and pipe 

1 flow {U mean flow velocity) case, respectively. In contrast to the Rayleigh- 
^^ \ Benard case the scaling exponents cannot (yet) be derived from the dynamical 

^ \ equations. 

^ ■ Pacs: 47.27.-i 

'NT ■ 

o^: 

O . I. INTRODUCTION 

o ■ 

■^ \ The relation between global flow properties and driving forces is interesting from a fun- 

• ^ ' damental point of view and for upscaling from laboratory experiments to applications. Ex- 

^ ■ amples are the change in mean flow through a pipe as a function of pressure drop, the 

>• . dependence of the heat transport as a function of temperature difference (Rayleigh-Benard 

KA \ (R-B) flow), and the increase in torque required to maintain a certain rotation speed in a 

5^ \ Taylor-Couette (TC) system. In all three systems the effects of the boundary layers are 

of prime importance but the way in which they are dealt with differs. For pipe flow the 

boundary effects are usually discussed in terms of the Prandtl-van Karman theory [jl| which 

assumes a logarithmic law for the proflle. [| For the connection between mean flow and 

pressure drop it predicts an implicit logarithmic relationship, the so called skin friction law 

)jl|], in reasonable agreement with the experimental data. 

On the other hand, Rayleigh-Benard (RB) convection has mainly been discussed in 
terms of algebraic relations: In the last decade the power law relation Nu ~ Ra?^"^ between 
the Nusselt number Nu and the Rayleigh number Ra was thought to be an appropriate 



c^ 



^Barenblatt and coworkers however also succeeded to describe the data in terms of power laws 



description of the experimental data [HH. Recently, it has turned out that the dependences 
are more involved 0-0- On the theoretical side, the analysis by Grossmann and Lohse 
Pi of the different dominant dissipation mechanisms leads to a detailed phase diagram 
for Rayleigh-Benard convection that is in good agreement with the latest 0,|^ and older 
experiments. In this theory the relations between Nu and Ra are again algebraic. 

For the Taylor-Couette system with rotating inner cylinder and resting outer one a 
description of the relation between the dimensionless torque G and the Reynolds number 
Re both in terms of a skin friction law and a pure power law has been tried ||^, ID]. The 



dimensionless torque G = T/pu'^L and the Reynolds number Re = Qa{b — a)/v are defined 
with T the torque, p the fluid density and v its kinematic viscosity, L the length of the 
cylinders, VL the angular rotation rate, and h and a the radii of the outer and inner cylinder, 
respectively. Lewis and Swinney's analysis [^ of their experimental data clearly shows that 
a pure power law G ~ -Re" with a = 5/3 as suggested in Refs. [|ll],0 does not describe the 



data. A description in terms of the skin friction law is in better agreement with the data. 
However, it still is not fully satisfactory, either, as the systematic drifts in their Fig. 1 show. 

It is our aim here to adopt Grossmann and Lohse 's Rayleigh-Benard theory to Taylor- 
Couette and pipe flow. In contrast to the RB case, it is presently not possible to derive the 
scaling exponents fully from dynamical equations. Instead, there will be one exponent that 
has to be and can be fitted consistently to data for both systems. 

We begin in section 2 with a discussion of the theory and the comparison to Lewis and 
Swinney's experimental data for TC flow. In section 3 we will adopt it to pipe flow for 



which precise high Re data on the pressure drop were obtained by Smits's group |]T3[. In 
section 4 we compare the quality of the data fit of this theory with that of the standard skin 
friction law theory. Section 5 gives conclusions. 

II. TAYLOR-COUETTE FLOW 

The basic idea ||^ behind the analysis of the thermal convection experiments is the 
splitting of the energy dissipation e into contributions from the boundary layer (BL) and 
the bulk, 

e = ^BL + ^bulk- (1) 

For the energy dissipation e in TC flow one strictly has 

u'^GQ 



27r(62 - a2) 



(2) 



This can be derived as usual by considering the energy balance in a Navier-Stokes flow. In 
analogy to ref. ^ (see there for an extensive discussion) etuik is estimated as 

(^buik ~ , "" ■ (3) 

b — a 

Here, f/^ is the typical velocity difference between the turbulent and the laminar (linear) 
profile. It is a measure for the turbulent activity induced by the rotation of the inner cylinder 
and it defines a Reynolds number Re^ = U^ip — a) /p. In the laminar case U^ = Rew = 0, 
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FIG. 1. a) Compensated plot G/i?e^/2+5?/2 ^g ^^1/2+^/2 ^^^^ ^ ^ -0.051. The points are 
Lewis and Swinney's data Q, the line the fit (^). Re varies in the range 10^ through 9 • 10^. 
b) Relative error (G — Gfit)/G of the suggested combination of power laws (|8|) with S, = —0.051 
(open boxes) and relative error (/ — ffit)/f of the friction law fit ( [Tl| ) (filled circles). 



and therefore ebuik = 0. Obviously, f/^ must not be confused with the velocity U = 2'KaQ of 
the inner cylinder or the corresponding Reynolds number Re = Qa{b — a) /v. The Reynolds 
number Re (or U) is imposed on the flow whereas Re^, (or Uui) is the response of the system. 
The situation can be compared with RB convection where the Rayleigh number is imposed 
on the cell whereas the response of the system is the large scale wind of turbulence, which 
again defines a wind Reynolds number Re^. Therefore, in analogy, also here we call U^ the 
wind velocity. 

It is this wind velocity which leads to the formation of a boundary layer of thickness A„. 
As in ref. [Bl we assume the BL to be of Blasius type [||], 



Am ~ 



{b-a)/^Re^. (4) 

For very small Re^ the BL will of course not diverge but saturate at a scale A^ ~ (& — a) 



which introduces different scaling relations [0 , but in the present work we are not interested 
in this very low Reynolds number regime. 

With its thickness A^ as the relevant length scale we estimate the energy dissipation in 
the BL as B 



eBL ~ ^ 



U,„ Xu 



Xlh 



Putting eqs. (0) to (|) together one obtains 

GRe = ciRel^^ + CiRe 



3 



(5) 



(6) 



where C\ and c^ are two unknown constants. The first term is the BL contribution, the 
second one the bulk contribution. 

The central question now is: How does Rtw depend on i?e? We do not know, but for 
large enough Re it seems reasonable to assume a power law dependence. 



Be 
~Re 



U 



Therefore, 



G = cii?e3/2+^«/2 + C2i?e2+3«. 



(7) 



(8) 



We perform a nonlinear fit of Lewis and Swinney's data |p to eq. (^, obtaining ^ = —0.051, 
ci = 10.5, and C2 = 0.196. The best way to check the quality of the fit (||) is to plot 
G jRe'l'^^^^l'^ vs Re^l"^^^!"^ so that according to eq. ^ a straight line should result. This 
indeed is the case, as shown in Fig. |l]a. The quality of the fit is underlined by the relative 
error shown in Fig. |I]b. 
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FIG. 2. a) Compensated plot y/i?e^/2+55/2 ^g ^gi/2+5/2 ^-^j^ ^ ^ -0.041. The points are data 
from ref. [13|, the line the fit (|To[). Re here varies between 3.16 • 10 and 3.53 • 10 . 
b) Relative error (p' — p'nt)lv' of the suggested combination of power laws ([lO|) with ^ = —0.041 
(open boxes) and relative error (/ — ffu)// of the friction law fit ( [ll| ) (filled circles). 



III. PIPE FLOW 



Consider now pressure driven flow through a pipe of radius R. The pressure gradient 
and the energy dissipation e per volume are related by 



Ap_ 
P 

4 



-u^ 



(9) 



where u^ =: U is the average of the x-velocity over the cross section of the pipe, which 
defines the Reynolds number Re = 2RU /v. As for the TC flow we define the wind velocity 
Uw as the maximal difference between the turbulent mean velocity profile and the laminar 
(parabolic) one. This maximum will occur close to the walls. Uw again defines a wind 
Reynolds number Rcu, = 2_Rf/^/i/. 

We split the energy dissipation in a boundary layer part and a bulk part as in equation 
(H) and estimate both contributions as above, eqs. @ and (||), with b — a replaced by R, and 
the thickness of the Blasius boundary layer being A„ ~ R/y/Re^. With this the equation 
for the dimensionless pressure drop p' = ApR^ / {pu'^) becomes 

p' = c;/?e3/2+5C/2 + c',Re^+^^, (10) 

where c[ and Cg are two unknown constants and ^ the power law exponent of eq. (^. 
Equation (|10D is the analog of equation (8) in the TC case. 



In order to test equation (|T0|) we consider the high precision pressure drop data of 
Smits and Zagarola [0. A nonlinear fit then results in ^ = —0.041 and c'l = 0.226 and 
4 = 0.00373. In Fig. |a we show p' / Re^/^+^^/^ vs i^e^/^+^Z^. If eq. (0) holds, a straight 
line should result which is the case. Again the quality of the fit is demonstrated with the 
relative errors in Fig. ^d. 

IV. THE SKIN FRICTION LAW 

Finally we would like to compare our description of the data with the standard skin 
friction law [|T|. Define the friction coefficient / = G/Re^ for TC flow and / = p'/Re^ for 
pipe flow. One then has, in both cases |ll|, |9|, p!0| , 

^ = c'MRe^-) + 4 (11) 

with two flow dependent constants c'/ and C2 which can be connected to the von Karman 
constant [|]. Employing equation (|Tl]) one can indeed fit both data sets reasonably. We do 
not show the fits as they have already been shown elsewhere (see figure 4a of ref. [^ for the 
TC case), but we present the relative error (/ — //jt)// in figures |l|a and |^ and compare it 
with the relative errors of the combined power laws eqs. (H) and (^^), respectively. For both 
TC and pipe flow the relative error of the friciton law is roughly of the same order as that 
of the fits (I) and (|10D , respectively, maybe somewhat larger. 

V. CONCLUSIONS 

The preceeding analysis shows that the splitting of the dissipation into a bulk and a 
boundary layer contribution as used in the Rayleigh-Benard theory can also be used to 



describe the Taylor-Couette flow and the pipe flow |]T3[ data. However, in contrast to the 
RB case, only one global balance equation is available, the one for the energy dissipation. 
Therefore, it is not possible to derive the asymptotic scaling exponents for both the wind 
Reynolds number and the dimensionless torque (dimensionless pressure drop) in the TC 



case (pipe case). Instead, one scaling exponent {^) must be fitted to the data. The ratio 
Uw/U ~ Re^ scales similarly in both cases, ^ = —0.051 and ^ = —0.041 for TC and pipe 
fiow, respectively. The relative error in both cases is less than one percent and within 
this precession the data indicate really different exponents. The origin of this difference is 
unclear. 

For shear fiow, the strict upper bound for energy dissipation is q = tL/U^ < 0.01087 for 
Re -^ oo [|l^,[l^. All the Couette experiments clearly lie below this bound, and even show 



a trend towards a scaling that is slower than U . When assuming the Kolmogorov length 



scale as smallest length scale on which dissipation contributes, Nicodemus et al. IT^ could 
numerically show that around Re = 10^ — 10^ one has c^ ~ i?e~°°^. Kerswell assumes the 
same cutoff and a certain background flow proflle and flnds c^ ~ Re~^''^ for Re — » oo [p!8 



Interestingly enough, the exponents 3^ = —0.12 (TC) and 3^ = —0.15 (pipe) that we flnd 
from the experimental data are very close to that value. However, the scaling Uw/U ~ Re^ 
with negative ^ implies that the wind does not increase as rapidly as the external velocity 
with Re and that according to our deflnition of U^ the relative difference between the laminar 
and the mean turbulent velocity proflles vanishes. Given the smallness of ^ the Reynolds 
numbers at which this could become signiflcant are not experimentally accessible. But the 
situation remains unsatisfactory and it clearly would be highly desirable to calculate this 
exponent more rigorously from the Navier-Stokes equations and to understand the relation 
to the mean flow proflle better. 
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